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Abstract: Drought was a severe recurring phenomenon in Iraq over the past two decades due to climate 
change despite the fact that Iraq has been one of the most water-rich countries in the Middle East in the 
past. The Iraqi Kurdistan Region (IKR) is located in the north of Iraq, which has also suffered from 
extreme drought. In this study, the drought severity status in Sulaimaniyah Province, one of four provinces 
of the IKR, was investigated for the years from 1998 to 2017. Thus, Landsat time series dataset, including 
40 images, were downloaded and used in this study. The Normalized Difference Vegetation Index (NDVI) 
and the Normalized Difference Water Index (NDWI) were utilized as spectral-based drought indices and 
the Standardized Precipitation Index (SPI) was employed as a meteorological-based drought index, to 
assess the drought severity and analyse the changes of vegetative cover and water bodies. The study area 


experienced precipitation deficiency and severe drought in 1999, 2000, 2008, 2009, and 2012. Study 


findings also revealed a drop in the vegetative cover by 33.3% in the year 2000. Furthermore, the most 
significant shrinkage in water bodies was observed in the Lake Darbandikhan (LDK), which lost 40.5% of 
its total surface area in 2009. The statistical analyses revealed that precipitation was significantly positively 
correlated with the SPI and the surface area of the LDK (correlation coefficients of 0.92 and 0.72, 
respectively. The relationship between SPI and NDVI-based vegetation cover was positive but not 
significant. Low precipitation did not always correspond to vegetative drought; the delay of the effect of 
precipitation on NDVI was one year. 
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1 Introduction 


Drought is one of the most complex natural disasters; it is still challenging to identify its cause 
and range of development (Aadhar and Mishra, 2017). In general, drought gives an impression of 
water scarcity due to insufficient precipitation, high evapotranspiration, and overexploitation of 
water resources, or a combination of all the above (Sheffield and Wood, 2008). Drought is 
triggered by an extended period of precipitation deficiency (Fadhil, 2011) and is usually 
associated with increased evaporation rates (Perez et al., 2016). Depending on how long and 
severe the lack of precipitation is, drought can bring various degrees of negative influence on 
soil/vegetation moisture, groundwater, rivers, streams, ecosystems, and human activities 
(Smakhtin and Hughes, 2004). 

Due to global warming, precipitation has decreased in many regions in the world (Yu et al., 
2019). Hence, this climate change is likely to cause more frequent and more severe droughts in 
several countries (Al-Quraishi et al., 2020). Various kinds of drought, such as agricultural, 
hydrological, meteorological, ecological, and socioeconomic drought, are caused by a reduction in 
the amount of water resources in areas characterized by a reduction in the availability of water 
resources (Sharma, 2006). Drought has severe impacts on agricultural and livestock production 
resources in any country. It negatively affects irrigated farmlands and areas that depend on rainfed 
irrigation farmlands. It also has serious consequences on population distribution, production, and 
life (Al-Quraishi et al., 2019; Gaznayee and Al-Quraishi, 20192). 

Geographically, Iraq is situated in the driest zone in the world; the water resources in Iraq are 
limited and mostly shared with other neighbouring countries (UNESCO, 2014). Water is the 
decisive factor for the stability and continuity of the agricultural department and the cornerstone 
of economic and social development in Iraq (Yaseen et al., 2018), which is heavily dependent on 
surface water and groundwater resources (Alobaidy et al., 2010). Unfortunately, Iraq has suffered 
from recurring droughts for the past 20 years (Street, 2012). 

Remote sensing (RS) technology can assess, monitor, and predict drought (Sun et al., 2010; 
Hazaymeh and Hassan, 2017). Proper indices derived from optical RS data have been used to 
assess and analyze drought severity (Hazaymeh and Hassan, 2017; Ghebrezgabher et al., 2019). 
RS and Geographic Information System (GIS) play vital roles in drought investigation as they 
offer practical and economical means to study the spatiotemporal distribution and evolution of 
vegetative cover and water resources through specific indices, such as the Normalized Difference 
Vegetation Index (NDVI) and the Normalized Difference Water Index (NDWIT) (Acharya and Ray, 
2007; Fadhil, 2013; Zhang et al., 2013). Based on the RS dataset, it is possible to employ spectral 
water indices to investigate the impacts of drought on water resources (Li et al., 2014; Gaznayee 
and Al-Quraishi, 2019b). The type of spectral indices for drought as well as their use depend on 
the kind of drought and the factors that lead to drought (Yang et al., 2014; Awchi and Jasim, 
2017). The Standardized Precipitation Index (SPI) is a well-known meteorological parameter 
widely used in drought studies (McKee et al., 1993; Edossa et al., 2010; Al-Quraishi and Negm, 
2020). Most drought studies often rely on meteorological and spectral indices so that appropriate 
measures and treatments can be established to mitigate the effects of drought (Loucks and van 
Beek, 2017; Mustafa, 2020). 

Although drought assessment using RS data and GIS technique is not a novel topic, to the best 
of our knowledge, few studies have employed these tools in the Iraqi Kurdistan Region (IKR) of 
Iraq, particularly in Sulaimaniyah Province. Therefore, we conducted this study for this reason. 
RS data were collected for the 20-year period from 1998 to 2017 (as this was a crucial period for 
the recurrence of drought episodes in Iraq) to calculate the NDVI and NDWI. The results were 
then integrated with meteorological data from the SPI and hydrological data (the surface area of 
the Lake Darbandikhan (LDK)) to analyze and assess drought trends in Sulaimaniyah Province in 
the IKR. 
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2 Materials and methods 


2. Study area 


This study was conducted in Sulaimaniyah Province in the IKR (Fig. 1), an area located in the 
north of Iraq. The IKR includes four provinces, with Sulaimaniyah Province being the largest one 
(two-thirds of the IKR). The study area is bordered by Iran from the east and surrounded by 
Kirkuk Province to the west, Salahalddin Province to the southwest, and Diyala Province to the 
south. Sulaimaniyah Province extends over 2.43x10^ km? and consists of 15 districts: 
Sulaimaniyah, Qaradagh, Sharazure, Saidsadiq, Penjwin, Halabja, Darbandikhan, Kalar, 
Khanagin, Kifri, Chamchamal, Dukan, Sharbazher (Mawat), Ranya, and Pishdar (Fig. 1a). LDK, 
one of the three largest lakes in the IKR, is situated 230 km northeast of Baghdad and nearly 60 
km southeast of Sulaimaniyah Province. LDK is connected to the Diyala River, a tributary of the 
Tigris River; it is a source of drinking and irrigation water for Sulaimaniyah Province and the 
surrounding areas, and has significant tourism benefits. 
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Fig. 1 Digital Elevation Model (DEM) of Sulaimaniyah Province (a) and locations of meteorological stations 
and geographical distribution of annual precipitation in Sulaimaniyah Province from 1998 to 2017 (b). District 
numbers: 1, Sulaimaniyah; 2, Qaradagh; 3, Sharazure; 4, Saidsadiq; 5, Penjwin; 6, Halabja; 7, Darbandikhan; 8, 
Kalar; 9, Khanaqin; 10, Kifri; 11, Chamchamal; 12, Dukan; 13, Sharbazher (Mawat); 14, Ranya; 15, Pishdar. The 
abbreviations of the districts are the same in Figure 3. 


2.2 Data collection 


Forty images of Landsat 5 TM, Landsat 7 ETM+, and Landsat 8 OLI representing the data (168/35 
and 168/36) for 20 years from 1998 through 2017 (Table S1) were assembled and used. The study 
period from 1998 to 2017 was chosen because it contained recurrent drought episodes. A mosaic of 
two Landsat scenes was constructed for each of the 20 years. The time series imageries were 
downloaded from the United States Geological Survey (USGS; https://glovis.usgs.gov/). All Landsat 
images were acquired in April or May in each of the 20 years, during which most of the vegetative 
proliferation occurred. Landsat images have a spatial resolution of 30 m. The Advanced Spaceborne 
Thermal Emission and Reflection Radiometer-Global Digital Elevation Model (ASTER-GDEM) V2 
dataset with a spatial resolution of 30 m, available from the National Aeronautics and Space 
Administration (NASA; https://www.nasa.gov/), was utilized as the Digital Elevation Model (DEM) 
for this study (Fig. 1a); the DEM was used to identify the relationships between vegetation 
percentage and altitude. 

Daily and monthly precipitation datasets at the ten meteorological stations in Sulaimaniyah 
Province were obtained from the Meteorological Department of the Ministry of Agriculture and 
Water Resources, Kurdistan Region Government (KRG), Iraq, during the period from 1998 to 2017. 
The detailed information of precipitation at the ten meteorological stations during the period from 
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1998 to 2017 is shown in Table S2. 
2.3 RS data processing 


After May 2003, the scan line corrector (SLC) of the Landsat 7 ETM+ sensor failed, resulting in 
images with strip line gaps, called SLC-off data type images. Those images had strip line gaps on 
their sides, causing approximately 22.0% data loss. To correct the SLC-off images and perform 
gap-filling, we employed the ENVI 5.3 (Harris Geospatial Solutions, Inc., USA) and a gap-fill 
plugin; it should be noted that it is not available in the ENVI Code Library currently, but a copy can 
be found at https://docs.google.com/file/d/OB3e_wo8OTO47b3c4ZHNyVONmUkk/edit? pli=1. 

This add-on module of the ENVI software platform fills gaps in one scene with data from another 
Landsat scene. The downloaded Landsat images were corrected by calibrating the digital number 
into radiance values using the data in Landsat image metadata files. Then, the resultant images were 
converted into top of atmosphere reflectance values using ENVI 5.3. The images were then 
georeferenced to the Universal Transverse Mercator, Zone 38 North, with the World Geodetic 
System 84 datum. A mosaic of two Landsat satellite scenes covering the entire study area was 
conducted for each of the 20 years of the study period. The produced mosaic represents and includes 
the entire study area. The image-to-image registration was accomplished with a root mean square 
error (RMSE) of 0.5 pixels. 


2.4 Spectral drought indices 
2.4.1 NDVI 


The NDVI is the most used index for vegetation monitoring. It accounts for all the green vegetation 
and is based on the combination of red band and near infra-red (NIR) band wavelengths and can be 
computed by the well-known formula (Eq. 1) of Rouse et al. (1974). 

NDVI-(NIR-RED)/(NIR* RED), (1) 
where NDVI is the Normalized Difference Vegetation Index; NIR is the near infra-red band 
(841—876 nm); and RED is the red band (620—670 nm). 

Healthy vegetation has less reflection in the visible range of the electromagnetic spectrum (EMS) 
due to the absorption of chlorophyll and other pigments. However, it has a high reflectance in the 
NIR section of the EMS. The NDVI is a strong vegetation signal and is mostly utilized to 
differentiate vegetative areas from non-vegetative areas (Huang et al., 2014). Its digital number 
values range from —1 to 1. Specifically, values from —1 to O represent the non-vegetative features, 
such as bare surface, built-up area, and water bodies, while values from 0 to 1 represent vegetative 
cover features. 


2.4.2 NDWI 


The NDWI is a parameter used for monitoring water body changes based on green band and NIR 
band wavelengths. It can be computed by the formula (Eq. 2) of McFeeters (1996): 

NDWI=(GREEN-NIR)/(GREEN+NIR), (2) 
where NDWI is the Normalized Difference Water Index; and GREEN is the green band (525-600 
nm). 

A water body has strong absorbability and low radiation from the visible spectrum to the infra-red 
range. According to McFeeters (1996), the water body can be mapped using a threshold value to 
better distinguish surfaces without detectable water (NDWI values less than 0.3) from those with 
detectable water (NDWI values higher than or equal to 0.3). In the current study, the NDWI was 
used to determine the size of water body, such as the LDK, and the NDWI fraction images were 
classified using the ArcGIS package. 


2.5 Meteorological drought index 


The SPI was proposed by McKee et al. (1993) and has been increasingly used during the last 
two decades because of its substantial theoretical development and robustness and versatility in 
drought analyses (Tsakiris et al., 2007; Nalbantis and Tsakiris, 2009). It can be calculated using 
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the following formula: 

SPI-(Xij-Xiy/oi, (3) 
where SPI is the Standardized Precipitation Index; Xj is the seasonal precipitation for the it 
station and Dh observation (mm); X;is the long-term mean precipitation for the i'^ station (mm); 
and ci is the standard deviation for the i station (mm). The SPI calculation for any location is 
based on the long-term precipitation record for the desired period. This long-term record is 
fitted to a probability distribution, which is then transformed into a normal distribution so that 
the mean SPI for the location and desired period is zero. The equation is widely used for 
monitoring meteorological drought and investigating and assessing its severity. It is also utilized 
to test the relationship of SPI with precipitation and crop production. 

One of the main advantages of SPI is that it only requires precipitation data for calculation. 
This simplicity makes it suitable in areas that lack diverse data, and renders it a popular tool. 
The SPI guarantees independence from a geographic location, and it can be calculated using the 
precipitation averages in the same location (Rossi et al., 2007). It was developed to measure the 
amount of precipitation decline for different time scales, reflect the reduced precipitation on 
water resource availability. Typically, calculating SPI requires precipitation data over a long 
period (at least 20 years). The precipitation time series dataset is fitted to a gamma distribution, 
and then it is transformed into a normal distribution with the aid of an equal probability 
transformation (Guttman, 1999). The function of probability density defines the gamma 
distribution: 


= 1 oli 

x,&, p= x e^, 4 
g(x, a, B) BTO (4) 
where a and £ are the shape and scale parameters, respectively; x is the amount of precipitation 
at a meteorological station (mm); and T (a) is the gamma function. The SPI calculation includes 
a gamma probability distribution fitted to a given frequency distribution for a meteorological 
station. For each meteorological station, the a and f parameters of the gamma probability 
density function were estimated for 12 months per year and for each year of the 20-year period. 

The maximum likelihood solutions were used to optimally estimate the a and f parameters: 


aa EE (5) 

f- X la, (6) 

A- In(X) Zi (7) 
n 


where n is the number of observations; X is the mean precipitation (mm); and A is the sample 
arithmetic and geometric means given by under some conditions. 

The results from Equation 8 were used to calculate the cumulative probability of an observed 
precipitation event. Since the gamma function is undefined for x=0 and a precipitation 
distribution may contain zero, the cumulative probability (H(x)) is the extended gamma: 

H(x)=q+ (1-4)xG@), (B 
where H(x) is the cumulative probability; q is the probability of zero; and G(x) is the cumulative 
probability of the incomplete gamma function. If m is the number of zero in a precipitation time 
series, q can be estimated by m/n. The cumulative probability was then transformed to the 
standard normal random variable (Guttman, 1999). The Drought Indices Calculator (DrinC) 
software used for calculating the SPI was designed by Tigkas et al. (2013). Positive SPI values 
indicate the availability of wet conditions, and the precipitation average is higher than the 
general precipitation average. In contrast, negative SPI values indicate the prevalence of dry 
conditions, and the precipitation average is lower than the general precipitation average in the 
region (Kamali et al., 2017). The anomaly strength was classified after the standardization of 
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SPI values, as shown in Table 1. Further, Figure 2 represents the general methodology flowchart 
adopted in this study, demonstrating the fundamental process of drought trend analysis. 


Table 1  SPI-based drought severity classes (McKee et al., 1993) 


SPI Classification 
22.00 Extremely wet 
1.50-1.99 Severe wet 
1.00-1.49 Moderately wet 
—0.99-0.99 Near normal 
—1.49— —1.00 Moderate drought 
—1.99- —1.50 Severe drought 
<-2.00 Extreme drought 


Note: SPI, Standardized Precipitation Index. 


Remote sensing obtaining, 
satallite 
images from USGS 


Meteorological data 
(rainfall data) for ten stations 


Data customization 
(excel sheet 
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Radiometric correction 
(DN>radiance>reflectance) 
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Clipping Darbandikhan lake area 
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Lake area calculation 
and color mapping 
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Vegetation cover calculation and 
color mapping in term of drought 
using ArcGIS 10.4.1 


Fig. 2 Flowchart of the methodology adopted in this study. NDVI, Normalized Difference Vegetation Index; 
NDWI, Normalized Difference Water Index; SPI, Standardized Precipitation Index; DN, digital number; RED, red 
band; GREEN, green band; NIR, near infra-red band; USGS, United States Geological Survey. 


3 Results 


All the indices (NDVI, NDWI, and SPI) and the surface area of the LDK were computed on a 
one-year basis due to limited imagery availability for the study area. 


3.1 NDVI variation 


The results shown in Table 2 and Figure 3 indicated that the most critical years for drought as 
detected from vegetation growth were 2000, 2008, 2009, and 2013. During these years, the 
vegetation cover was remarkably reduced compared with that in other years. Drought was most 
severe in 2000, and the vegetative cover area reduced to 3.31x10? km? (occupying 13.6% of the total 
study area). The average of the vegetation coverage is 49.6% during 1998-2017, and the vegetation 
coverage in 2000 fluctuated 33.396 from the average. 

Figure 3 shows the spatiotemporal distribution of NDVI in Sulaimaniyah Province from 1998 to 
2017. The vegetation cover had a considerable spatial variability at the spatial scale, mainly in the 
center of Sulaimaniyah Province, whereas the northeastern and southern parts constantly remained 
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the most and least vegetative areas, respectively. For the study period from 1998 to 2017, there was 
a strong relationship between NDVI and precipitation averages across most parts of Sulaimaniyah 
Province. 


Table 2 Characteristics of NDVI and vegetative cover in the study area from 1998 to 2017 


NDVI Vegetative cover 
id Max. Min. Mean SD Area (x10? km?) Coverage (96) Fluctuation (96) 
1998 0.992 0.098 0.242 0.120 10.09 41.5 DA 
1999 0.983 0.098 0.230 0.101 12.17 50.1 3.1 
2000 0.553 0.016 0.097 0.098 3.31 13.6 —33.3 
2001 0.724 0.029 0.201 0.126 14.78 60.5 13.6 
2002 0.726 0.064 0.233 0.118 13.58 55.9 9.0 
2003 0.722 0.051 0.245 0.128 11.95 49.2 2.3 
2004 0.674 0.039 0.204 0.112 12.23 50.3 3.4 
2005 0.668 0.057 0.167 0.101 10.84 44.6 -2.3 
2006 0.782 0.018 0.214 0.143 14.99 61.7 14.8 
2007 0.730 0.108 0.318 0.127 12.93 53.2 6.3 
2008 0.611 0.019 0.116 0.081 8.88 36.5 -10.4 
2009 0.696 0.084 0.234 0.110 8.84 36.4 —10.5 
2010 0.655 0.127 0.258 0.097 10.71 44.1 -2.8 
2011 0.598 0.061 0.150 0.068 11.09 45.6 -1.3 
2012 0.691 0.013 0.211 0.129 11.85 48.8 1.8 
2013 0.628 0.163 0.267 0.079 9.37 38.6 -8.4 
2014 0.862 0.288 0.484 0.125 13.47 55.4 8.5 
2015 0.642 0.182 0.320 0.085 13.43 55.3 8.3 
2016 0.646 0.175 0.289 0.078 12.87 53.0 6.1 
2017 0.640 0.178 0.278 0.070 10.75 44.2 -2.7 


Note: NDVI, Normalized Difference Vegetation Index; Max., maximum; Min., minimum; SD, standard deviation. Fluctuation in the last 
column indicates fluctuation around the mean vegetation coverage. 


3.2 SPI variation 


Table S3 and Figure 4 show the calculated SPI values for Sulaimaniyah Province from 1998 to 
2017. The results showed an irregular cyclical pattern of dry/wet periods during the study period. 
Droughts were observed in hydrological years of 1998-1999, 1999-2000, 2007-2008, 
2008-2009, and 2011-2012, with some differences across the meteorological stations. Generally, 
droughts commonly occurred at the beginning of the rainy season, reflecting either a delay in 
rainfall or a decrease in rainfall. However, the highest level of drought was observed during 
1998-1999 in Sulaimaniyah, Bazian, Halabja, Penjwen, Saidsadiq, and  Darbandikhan 
meteorological stations, with average SPI values around —2.22, —1.90, —1.94, —2.08, —1.48, and 
—1.72, respectively (Table S3; Fig. 4). Moreover, Bazian, Halabja, Darbandikhan, Chamchamal, 
Kalar, and Pebaz meteorological stations also experienced moderate to severe droughts during 
2007-2008, with average SPI values of —2.38, —1.92, —1.92, —2.17, —2.56, and —2.08, respectively 
(Table S3; Fig. 4). 

On the other hand, the highest degree of drought severity occurred in Kalar, where the average 
SPI value was —2.56 during 2007—2008 (Table S3; Fig. 4). Figure 5 shows that SPI fluctuated in 
Darbandikhan meteorological station in Sulaimaniyah Province from 1998 to 2017. The first SPI 
value decline began in 1999 and continued until 2001. This reduction (—1.72) is consistent with 
the decrease in precipitation observed in Sulaimaniyah Province throughout 1998-1999, as it was 
a dry year (Table S3; Fig.4). The second SPI value decline began in 2007 with the decrease in 
precipitation in Sulaimaniyah Province throughout 2007-2008 (Table S3; Fig. 4). 
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Fig.3 Spatiotemporal variations of NDVI-based vegetation cover in the study area from 1998 to 2007 
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Fig. 5 Annual precipitation, SPI, and the surface area of the LDK (Lake Darbandikhan) at Darbandikhan 
meteorological station from 1998 to 2017 
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3.3 LDK area and NDWI variation 


The combination of NDWI and SPI is a very effective method for investigating drought patterns to 
monitor the surface area of the LDK and the river water levels. We found that the effect of drought 
on the surface area of the LDK was severely pronounced in 1999, 2000, 2001, 2008, and 2009, 
based on NDWI values (Table 3; Fig. 6). The mentioned years had lower SPI values, less surface 
area of the LDK, and lower precipitation averages. 

The surface area of the LDK changed throughout the studied 20 years (Table 3; Fig. 6). 
Specifically, the largest surface area was observed in 2003, 2005, and 2016, with the values of 
113.70, 113.50, and 114.40 km?, respectively; the lowest surface area was found in 1999 and 2009, 
with the values of 37.50 and 39.00 km», respectively. 


Table3 Surface area of the LDK and its percentage change from 1998 to 2017 


Year Surface area (km?) Fluctuation (96) 
1998 103.60 25.6 
1999 37.50 —40.5 
2000 49.50 —28.5 
2001 48.60 -29.4 
2002 93.80 15.8 
2003 113.70 35.7 
2004 94.50 16.5 
2005 113.50 35.5 
2006 93.10 15.1 
2007 89.00 11.0 
2008 49.72 —28.3 
2009 39.00 —39.0 
2010 66.40 -11.6 
2011 53.40 —24.6 
2012 92.30 14.3 
2013 81.40 3.4 
2014 60.10 -17.9 
2015 59.00 —19.0 
2016 114.40 36.4 
2017 106.90 28.9 


Note: LDK, Lake Darbandikhan. Fluctuation indicates the fluctuation around the average percentage of the surface area of the LDK. 


Although the IKR has abundant water resources, these resources are generally considered to 
be limited and to change with time and place. Based on data from the Ministry of Agriculture 
and Water Resources in Iraq, we found that about 40.096 of the region's spring water dried up 
during previous drought events (UNESCO, 2014). 


3.4 Pearson correlation analysis 


The Pearson correlation analysis between NDVI and SPI showed that NDVI had a low and 
insignificant correlation coefficient with SPI (Table 4), which indicated that vegetation cover 
extension in a specific year (assessed from April to May) was not dependent on the rainfall (or 
SPI) measured in the same year. Further, there was a significant positive relationship (P«0.01) of 
SPI with precipitation (r=0.92) and surface area of the LDK (r=0.73). These results are consistent 
with those of Gaznayee and Al-Quraishi (2020). 


4 Discussion 


The lowest NDVI values were found in the southern and western regions of the study area; the 
northeastern region constantly had more developed vegetation and higher NDVI values (Fig. 3). 
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There was a significant decline in the annual precipitation averages in some districts in 
Sulaimaniyah Province in 2000, 2008, and 2009 (Table S2). The vegetation coverage percentage 
of the NDVI-based vegetative areas shown in Table 2 was due to the fluctuation in the 
precipitation averages shown in Table S2. Furthermore, low precipitation and high temperature 
played significant roles in decreasing NDVI values and vegetative areas in the southwestern parts 
of Sulaimaniyah Province during the growing seasons, which are depicted in the correlation 
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Fig.6 Spatial distribution changes of the LDK surface area from 1998 to 2017 
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Table A Pearson correlation coefficients among the NDVI, precipitation, SPI, and the surface area of the LDK 


NDVI Precipitation SPI Surface area of the LDK 
NDVI 1.00 0.20 0.24 0.13 
Precipitation 0.20 1.00 0.92" 0.72" 
SPI 0.24 0.92* 1.00 0.73" 
Surface area of the LDK 0.13 0.72" 0.73" 1.00 


Note: ™ indicates that correlation is significant at P«0.01 level. 


coefficient statistics. The consecutive spatiotemporal variations of NDVI can reveal the onset and 
extent of drought over the studied 20 years. Based on the results of NDVI variation shown in 
Table 2 and Figure 3, we conclude that Sulaimaniyah Province experienced drought episodes 
mainly in 2000, 2008, and 2009. The results mentioned above are consistent with Almamalachy 
(2020), who found that the year 2000 was characterized by the most severe drought compared to 
other years during 1998-2017. 

Results of SPI variation indicated that the severity of drought was profound during 
1998-1999 and 1999-2000; in 2007—2008, and the southeastern region was the most affected 
area. The western and southern parts of Sulaimaniyah Province were influenced by moderate 
and severe drought in 2012 (Table 1 and S3). The hydrological years of 2001—2002, 2009-2010, 
and 2015-2016 were wet compared with other hydrological years (Figs. 4—5). According to the 
study of McKee et al. (1993), drought occurs when the SPI is negative, while there is no drought 
occurring when the SPI is positive. 

The SPI values did not have a noticeable trend over the study period (Table S3). However, there 
were two periods when the drought was severe: period from 1998-1999 to 1999-2000 and period 
from 2007-2008 to 2008-2009 (Figs. 4—5). The precipitation deficits lasted for at least three 
years; hence, the drought that occurred during these two periods can be considered long-term 
drought. The SPI results indicated that low precipitation does not always correspond to vegetative 
drought (Figs. 3—4). A positive relationship existed between SPI and NDVI, especially in drought 
years of 2000, 2008, and 2013. However, vegetation cover extension in a specific year (assessed 
from April to May) was not dependent on the precipitation measured in the same year. This is 
likely due to the delayed effect of precipitation on vegetation vitality. Ji and Peters (2003) 
reported that the effect of rainfall on vegetation became detectable more than one month later. A 
similar study was done by Wang et al. (2003), who found that the delay of the response of NDVI 
to precipitation was 4—8 weeks. In this study, the delay of the effect of precipitation on NDVI was 
one year. A low SPI in 1998-1999 caused a low NDVI in 2000, a low SPI in 2007—2008 caused a 
low NDVI in 2009, and a low SPI in 2011-2012 caused a low NDVI in 2013 (only in the western 
and southern parts of the study area). The lack of coordination between SPI and hydrological and 
agricultural drought is more prominent during the growing seasons (Bhuiyan et al., 2006). 

The decrease in the surface area of the LDK is consistent with the shrinkage in the lake surface 
area reported by UNESCO (2014). Fadhil (2011) investigated the impacts of drought on the LDK 
in terms of climate change, and found that the LDK experienced the most significant drought 
events in 2008 and 2009. The decline in surface area of the LDK could be attributed to the 
reduction in the amount of water flowed into Iraq from neighboring countries, since 
approximately 70.096 of the water in the LDK originates from Iran. Furthermore, fluctuation in 
the amount of water from rivers of neighboring countries such as Turkey and Iran also caused 
increasing and frequent droughts in the IKR. In this study, based on the NDWI values, we found 
that the effects of severe droughts on the surface area of the LDK were pronounced in 1999 and 
2009, the years characterized by low precipitation and SPI values. 


5 Conclusions 


The present study investigated the occurrence and effects of drought in Sulaimaniyah Province in 
the IKR, over a period of 20 years (1998-2017). Landsat optical remote sensing images were 
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collected in the same season (April-May) during each year, and meteorological data were also 
used. Three standard spectral indices, i.e., SPI, NDVI, and NDWI, were evaluated on a one-year 
basis to assess the drought and its effect on vegetation cover and water resources with respect to 
the surface area of the LDK. Drought was more severe in 1999-2000 and 2008-2009, the years 
that were characterized by low SPI values and low surface area of the LDK due to reduced 
precipitation. Considering the significant recurrence of drought, it is crucial to satisfy the water 
needs of the study area by using other available water resources, such as snowfall, and 
groundwater, to increase the amount of water flow into the LDK. 

Statistical analyses highlighted that SPI is highly correlated with the surface area of the LDK. 
The correlation between SPI and NDVI in the same measured year was not significant, which was 
likely due to the delayed effect of scarce precipitation on vegetation. More detailed investigations 
are needed to understand the frequency of drought and its relationship to factors affecting itself. 
Additionally, efforts should be made to understand the effects of drought on crop production, 
forests, sustainability, water bodies, and socioeconomic resources. Due to the fact that only a 
one-year level analysis was conducted in this study, the jointly combined research of SPI and 
NDVI on a shorter-term basis (seasonal or even monthly level), in order to highlight the delay of 
the direct effect of drought on vegetation cover, and therefore on NDVI, is needed. We advise that 
Google Earth Engine may be a possible resource to develop such an analysis in the future. 
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Appendix 


Table S1 Information of Landsat 5 TM, Landsat 7 ETM+, and Landsat 8 OLI images in the study area from 
1998 to 2017 


Year Sensor Path/Row Acquisition date (dd/mm) Resolution (m) 
1998 Landsat 5 TM 168/35, 168/36 30/05, 30/05 30 
1999 Landsat 5 TM 168/35, 168/36 01/05, 01/05 30 
2000 Landsat 7 ETM+ 168/35, 168/36 25/04, 25/04 30 
2001 Landsat 7 ETM+ 168/35, 168/36 28/04, 28/04 30 
2002 Landsat 7 ETM+ 168/35, 168/36 01/05, 01/05 30 
2003 Landsat 7 ETM+ 168/35, 168/36 20/05, 20/05 30 
2004 Landsat 7 ETM+ 168/35, 168/36 06/05, 06/05 30 
2005 Landsat 7 ETM+ 168/35, 168/36 23/04, 23/04 30 
2006 Landsat 7 ETM+ 168/35, 168/36 12/05, 28/05 30 
2007 Landsat 7 ETM+ 168/35, 168/36 07/05, 07/05 30 
2008 Landsat 7 ETM+ 168/35, 168/36 15/04, 15/04 30 
2009 Landsat 7 ETM+ 168/35, 168/36 20/05, 20/05 30 
2010 Landsat 7 ETM+ 168/35, 168/36 05/04, 19/04 30 
2011 Landsat 5 TM 168/35, 168/36 16/04, 15/04 30 
2012 Landsat 7 ETM+ 168/35, 168/36 26/04, 26/04 30 
2013 Landsat 8 OLI 168/35, 168/36 23/05, 23/05 30 
2014 Landsat 8 OLI 168/35, 168/36 24/04, 24/04 30 
2015 Landsat 8 OLI 168/35, 168/36 27/04, 27/04 30 
2016 Landsat 8 OLI 168/35, 168/36 15/05, 15/05 30 


2017 Landsat 8 OLI 168/35, 168/36 18/05, 18/05 30 
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